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OBJECTIVE: 

The  electrodynamic  mechanisms  of 
blackbody  emission  of  infrared  (IR) 
radiation  from  matter  are  well  known,  as 
testified  by  the  prolific  literature  existing  in 
the  field  and  covering  scientific  aspects  and 
applications  [1-3].  The  reversed  process, 
however,  i.e.  the  interaction  of  IR  radiation 
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with  matter,  is  still  obscure.  This  Project  aims  at 
unveiling  the  basic  mechanisms  involved  in  the  IR 
radiation-matter  interaction,  and  outlining  the 
possible  applications  related  to  harvesting  IR 
radiation  as  a  sustainable  and  renewable  source  of 
energy.  These  two  goals  are  achieved  by 
illuminating  a  power  generator  (PG)  device, 
shown  in  Fig.  1,  with  IR  radiation,  and  studying 
the  voltage  difference  AV,  the  current  I,  the  power 
AV*I,  and  temperature  difference  AT  produced  as 
a  function  of  time  t  upon  starting  the  illumination. 

This  procedure  is  IR  power  generation.  In 
addition,  this  Project  also  aims  at  (i)  training 
undergraduate  students  in  energy-related  science 
and  technology  relevant  for  the  U.S.  Navy,  (ii) 
disseminating  knowledge  derived  from  IR  power 
generation  to  general  and  scientific  audience,  and 
(iii)  developing  knowledge  relevant  for  a  course  in 
Thermodynamics  and  Statistical  Mechanics. 

Scientific  and  technological  objectives.  In 
Year-1  we  discovered  the  existence  of  an  electric  and  an  cntropic  contribution  in  the  energy 
transfer  from  IR  radiation  to  matter  through  the  PG  device  [4],  Additionally,  we  found  that,  with 
the  electric  contribution,  the  PG  device  has  a  capacitor-type  behavior,  whereas  it  has  a 
thermoelectric-type  behavior  with  the  entropic  contribution.  We  also  started  examining  the 
effects  on  the  voltage  difference  AV(t)  of  (i)  colored  tape  covering  the  illuminated  face  of  the  PG 
device  (Fig.  2),  (ii)  angle  of  incidence  0  of  the  IR  radiation  with  respect  to  the  PG  device  (Fig. 
3),  (iii)  in-series  coupling  of  two  PG  devices  (Fig.  3),  (iv)  variable  laser  power  [with  Marc  Currie 
at  the  Naval  Research  Laboratory  (NRL)  and 
Oleksandr  Kokhan  (JMU)]  (Fig.  4),  (v)  nano-scale 
architecture  of  the  PG  device  (Figs.  5  and  6)  with 
active  element  consisting  of  a  thin  film  of  variable 
resistance  fabricated  via  atomic  layer  deposition 
(ALD)  by  Virginia  D.  Wheeler  at  NRL.  We 
disseminated  the  knowledge  we  acquired  through  the 
Workshop  on  “ Infrared  radiation,  Thermoelectricity 
and  Chaos ”  which  took  place  on  June  17,  2015  at 
JMU.  The  Workshop  was  funded  by  the  ONR 
Award  #  N000141410378.  In  Year-2  we  will 
complete  the  objectives  of  Year-1.  In  addition,  we 
will  (i)  develop  a  macro-PG  device  with  Boris 
Feygelson  at  NRL,  (ii)  study  the  current  I(t) 
produced  when  the  PG  device  is  illuminated  by  IR 
radiation,  and  (iii)  model  the  behavior  in  time  t  of  the 
variables  at  play  in  IR  power  generation  adopting  the 
conjugated-variable  based  approach  to  capture  the 
mechanism  triggering  the  switch  in  behavior  of  the 
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Fig.  3.  (a)  Amplitude  of  jump  in  voltage, 
A(AV),  immediately  after  starting  the 
illumination  with  IR  radiation  of  a  PG 
device  in  various  settings  (tape  or  in¬ 
series  coupling)  versus  angle  of  incidence 
0.  (b)  A(AV)  obtained  when  a  PG  device 
is  activated  by  a  100  LI  resistor  placed  in 
contact  with  the  active  face  of  the  device 
in  various  device  geometries. 
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Fig.  2.  (a)  Amplitude  of  jump  in  voltage, 
A(AV),  immediately  after  starting  the 
illumination  with  IR  radiation  of  a  PG 
device  in  geometries  named  A  and  T  and 
covered  or  not  with  tape.  (b)  A(AV) 
obtained  when  a  PG  device  is  activated  by  a 
100  Q  resistor  placed  in  contact  with  the 
active  face  of  the  device  in  geometries  A 
and  T,  and  covered  or  not  with  tape  [4]. 
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PG  device  from  capacitor-type  to 
thermoelectric-type ,  and  vice-versa. 

Educational  objective.  In  Year-1,  (i)  7 
undergraduate  students  participated  in  research 
and  dissemination  activities  related  to  this 
Project,  and  (ii)  we  shared  knowledge  and 
awareness  on  IR  power  generation  to  a  general 
audience  [5],  and  the  broader  scientific 
community  [6-12].  In  Year-2  we  plan  on 
disseminating  our  results  in  venues  such  as 
AVS  63rd  International  Symposium,  94th 
Annual  Meeting  of  the  Virginia  Academy  of 
Science,  2016  Undergraduate  Research 
Symposium  in  Physics  (JMU),  2016  Summer 
REU  Research  Symposium  (JMU),  and, 
possibly,  SPIE  Optics  and  Photonics 
Nanoscience  and  Engineering  Conference;  San 
Diego  (CA),  August  30-September  1, 2016.  In 
addition,  we  plan  on  developing  an  innovative  approach  to  Thermodynamics  and  Statistical 
Mechanics  using  a  conjugated-variahle  based  approach.  We  showed  that  this  concept  is  useful 
in  recognizing  the  electric  and  the  entropic  contributions  in  IR  power  generation  [4],  and  propose 
that  it  might  be  useful  also  for  other  physical  phenomena.  We  plan  on  drafting  a  textbook  in 
Thermodynamics  and  Statistical  Mechanics  possibly  also  with  additional  funding. 

Significance  and_  Future  Naval  Relevance;.  The  low  frequency  and  low  power  radiation 
studied  in  this  Project  is  important  for  fundamental  studies  and  also  for  applications  on  ships, 
satellites,  cars,  personal  backpacks,  and,  more  generally,  where  non-dangerous  energy  is  needed 
at  all  hours  of  the  day,  independent  of  weather  conditions.  Full  understanding  of  power  AV*I 
generation,  and  designing  robust  devices  (nano-PG  or  macro-PG  devices)  are  important 
objectives  to  implement  the  concept  of  IR  power  generation  and  move  it  into  applications.  The 
importance  of  our  studies  lies  also  in  that  the  method  of  harvesting  IR  radiation  using  a  PG 
device  could  be  extended  to  other  electromagnetic  radiation.  Broadening  the  spectrum  of  the 
usable  electromagnetic  radiation  would  greatly  contribute  to  the  exploitation  of  the  electric 
contribution  for  an  efficient  sustainable  and  renewable  source  of  energy  available  to  humankind. 

TECHNICAL  REPORT. 

Scientific  achievements. 

(i)  Decoupling  of  the  electric  and  entropic  contributions  in  IR  power  generation.  We 

discovered  that  energy  transfers  from  IR  radiation  into  matter  (e.g.  the  illuminated  surface  of  the 
PG  device)  affects  the  couples  of  conjugated  variables  consisting  of  charge  q  and  voltage  V,  on 
one  hand,  and  entropy  E  and  temperature  T,  on  the  other.  The  q-V  couple  generates  the  electric 
contribution  and  dominates  when  IR  radiation  activates  the  PG  device  [4].  This  contribution 
imparts  on  the  PG  device  a  capacitor-type  behavior,  and  is  sensitive  on  how  the  insulating  and 
conducting  layers  alternate  in  the  PG  device  (Fig.  1).  The  E-T  couple  generates  the  entropic 
contribution,  and  dominates  when  Joule  heating,  due  to  an  electrical  resistance,  activates  the  PG 
device  [4],  This  contribution  imparts  on  the  PG  device  a  thermoelectric-type  behavior  obeying 
the  Seebeck  effect:  AV(t)  =  S*AT(t),  where  S  is  the  Seebeck  coefficient.  We  have  yet  to  explain 


Average'Power  [rriW] 


Fig.  4  Amplitude  of  jump  in  voltage,  A(AV), 
and  temperature,  A(AT),  immediately  after 
starting  the  illumination  with  an  IR  laser  of  a 
PG  device  versus  average  power  for  (a)  CW 
laser,  and  (b)  pulsed  laser.  Panel  (c)  zooms  into 
the  low  power  region  of  panel  (b). 
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what  switches  the  PG  device’s  behavior  from  capacitor- 
type  to  thermoelectric-type ,  and  vice-versa,  and  plan  on 
using  the  conjugated-variable  based  approach  in  this 
effort.  Relevant  couples  of  conjugated  variables  will  be 
charge  q  and  voltage  V,  and  entropy  E  and  temperature 
T.  Here,  q=G(z,t)*surface  area,  where  G(z,t)  is  surface 
charge  density. 

(ii)  Effect  of  colored  tape  on  the  production  of 
voltage  AV(t)  in  IR  power  generation.  We  discovered 
that  when  the  electric  contribution  is  at  play  in  the 
interaction  between  IR  radiation  and  matter  in  a  PG 
device,  the  amplitude  of  the  voltage  jump  A(AV) 
achieved  upon  starting  the  illumination  can  be  increased 
at  least  twice  by  simply  placing  tape  on  the  face  of  the 
PG  device  hit  by  IR  radiation.  The  finding  is  illustrated 
in  Fig.  2a  [4],  and  is  not  detected  when  the  entropic 
contribution  dominates,  as  can  be  seen  in  Fig.  2b.  This  result  (a)  decouples  the  electric  from  the 
entropic  contribution,  and  (b)  suggests  that  tape  can  modify  the  capacitance  C  of  the  PG  device. 
Thus,  we  need  to  investigate  on  the  role  of  C  in  switching  the  PG  device  behavior’s  from 


Fig.  5.  Schematic  of  a  nano-PG 
device  illuminated  by  IR  radiation. 
The  ALD  film  connects  electrically 
the  “hot”  and  “cold”  faces  of  the 
device  made  of  copper  plates. 


A(AV)  (Tape}/A(AV)  =2.5 


capacitor-type  to  thermoelectric-type,  and  vice-versa. 

(iii)Effect  of  angle  of  incidence  0  on  the  production  of  AV(t)  in  IR  power  generation. 

We  found  that  the  amplitude  of  the  voltage  jump  A(AV)  slightly  increases  with  the  angle  of 
incidence  0  of  the  IR  radiation  on  the  PG  device.  The  phenomenon  is  illustrated  in  Fig.  3a.  We 
did  not  observe  a  clear  dependence  on  0  of  A(AV)  when  the  PG  device  is  activated  by  Joule 
heating,  as  shown  in  Fig.  3b.  With  will  need  to  collect  more  data  to  better  characterize  this 
phenomenon,  and,  possibly,  use  it  to  decouple  the  electric  from  the  entropic  contribution  in  IR 


A(AV)  {45°}/A(  AV)  =1.2 


power  generation. 

(iv) Effect  of  in-series  coupling  of  two  PG  devices  on  the  production  of  AV(t)  in  IR 
power  generation.  We  coupled  two  PG  devices  in  series,  and  found  that  the  amplitude  of  the 
voltage  jump  A(AV(t))  increases  after  starting  the  illumination  with  IR  radiation,  as  can  be  seen 
in  Fig.  3a.  The  amplitude  can  increase  up  to  three  times  when  tape  is  used  in  addition  to  in¬ 
series  coupling.  We  did  not  observe  a  clear  increase  in  A(AV)  when  the  PG  devices  are  activated 
by  Joule  heating,  as  shown  in  Fig.  3b.  Also  in  this  case  we  will  need  to  collect  more  data  to 
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exploit  this  finding  and  use  it  to  decouple  the  electric 


generation. 
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Fig.  6.  Amplitude  of  jump  in  voltage, 
A(AV),  and  temperature,  A(AT), 
immediately  after  starting  the  activation 
of  a  nano-PG  device  with  IR  radiation 
and  a  100  O  resistor  placed  in  contact 
with  the  active  face  of  the  device. 


(v)  Effect  of  variable  laser  power  on  the 
production  of  AV(t)  in  IR  power  generation.  In 

collaboration  with  Marc  Currie  at  NRL  and  Oleksnadr 
Kokhan  at  JMU,  we  are  testing  the  effects  of  IR 
radiation’s  power  on  PG  devices.  We  are  using  monochromatic  lasers  at  wavelength  A,=1064  nm 
with  range  of  average  power  from  9  to  700  mW.  With  continuous  wave  (CW)  lasers,  the  power 
remains  constant  throughout  the  duration  of  the  illumination.  With  pulsed  lasers,  the  power 
changes  from  a  minimum  to  a  peak  value  over  a  well-defined  pulse  duration.  The  peak-power  in 
our  experiments  spans  from  the  kW  to  the  MW  range.  Our  results  with  the  CW  lasers, 
summarized  in  Fig.  4a,  indicate  that  the  amplitude  of  the  voltage  jump  A(AV)  and  of  the 
temperature  jump  A(AT)  increase  linearly,  however  the  slope  of  A(AV)  vs.  average  power  is 
0.013  mV/mW,  one  order  of  magn  tude  larger  than  that  of  A(AT)  vs.  average  power,  which  is 
0.003  °C/mW.  Wc  interpret  this  discrepancy  as  an  effect  of  the  prevalence  of  the  electric 
contribution,  and  of  its  independence  from  the  entropic  contribution.  On  the  other  hand,  our 
results  for  the  pulsed  laser  case,  summarized  in  Fig.  4b  and  4c,  indicate  that  the  amplitudes  of 
A(AV)  and  A(AT)  increase  linearly  vs.  average  power  and  with  similar  slopes,  0.008  mV/mW 
and  0.009  °C/mW,  respectively.  This  behavior  is  consistent  with  the  Seebeck  effect,  and  points 
out  the  prevalence  of  the  entropic  contribution.  From  the  experimental  point  of  view,  we  will 
need  to  investigate  :he  effects  of  (i)  colored  tape,  (ii)  angle  of  incidence  0,  and  (iii)  in-series 
coupling  of  PG  devices  illuminated  oy  the  lasers  to  gain  further  insight  on  what  causes  switch  of 
the  of  PG  devices  behavior  from  capacitor-type  to  thermoelectric-type  when  illuminated  by  CW 
and  pulsed  lasers,  respectively. 

(vi) Effect  of  a  nano-scale  architecture  of  the  PG  device  with  active  element  consisting  of 
a  thin  film  layer  fabricated  via  ALD  with  variable  resistance.  With  the  collaboration  of 
Virginia  D.  Wheeler  at  NRL  we  have  fabricated  thin  films  of  variable  resistance  to  build  nano- 
PG  devices  as  sketched  in  Fig.  5,  and  test  their  response  to  IR  radiation  and  Joule  heating.  The 
thin  films  so  far  consist  of  stacks  of  conducting  (e.g.  TiN)  and  non-conducting  (e.g.  TiCb)  layers 
deposited  via  ALD.  This  technique  enables  precise  thickness  and  doping  control.  The  thin  film 
resistance  so  far  spans  from  about  90  Q/sq  to  67  k£2/sq.  Our  results,  summarized  in  Fig.  6, 
indicate  that  the  illumination  of  the  nano-PG  device  with  IR  radiation  raises  the  amplitude  of  the 
voltage  jump  A(AV)  with  increasing  the  thin  film  resistance.  This  finding  is  quite  unexpected 
and  not  yet  understood.  We  plan  on  exploring  whether  the  observed  trend  persists  at  even  higher 
and  lower  resistances. 
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(•)  Conclusions:  The  current  picture  of  the  interaction  of  IR  radiation  with  matter,  achieved 
through  our  studies  on  IR  power  generation,  is  as  follows.  The  transfer  of  energy  from  IR 
radiation  obeys  the  laws  of  conservation  and  is  constrained  by  the  balance  between  the  following 
couples  of  conjugated  variables  active  in  the  process:  charge  (q)  and  voltage  (V),  and  entropy 
(E)  and  temperature  (T).  When  low  power  IR  radiation  illuminates  a  PG  device,  the  electric 
contribution  dominates  and  dictates  a  capacitor-type  behavior  to  the  device.  In  this  situation,  q 
and  V  are  related  by  the  capacitance  C  of  the  PG  device  and,  by  properly  engineering  C, 
maximum  efficiency  can  be  achieved  in  harvesting  IR  radiation.  When  the  power  of  the  IR 
radiation  is  increased  to  the  kW  range  and  above,  the  entropic  contribution,  related  to  E  and  T, 
dominates  and  dictates  a  thermoelectric-type  behavior  to  the  PG  device.  This  situation 
suppresses  voltage  production  at  the  advantage  of  temperature  increase.  The  circumstances  of 
the  switch  in  the  PG  device  behavior  must  be  investigated  and  might  involve  an  entanglement 
among  the  conjugated-variables  q-V  and  E-T.  These  findings  are  important  for  harvesting  IR 
radiation  as  a  renewable  and  sustainable  source  of  energy.  The  picture  on  the  behavior  of  the 
conjugated  variables  might  impact  other  scenarios:  e.g.,  in  biology  and  medicine  the  nerve 
system  is  viewed  as  capacitor  system  which  might  be  affected  by  low  power  IR  radiation  [13, 
14]  in  ways  similar  to  those  described  in  this  research. 

Educational  achievements. 

The  educational  results  achieved  are  as  follows: 

(•)  Student  training;.  In  Summer  2015  three  undergraduate  students  (Justin  M.  Kaczmar, 
Graham  P.  Gearhart  and  Tara  R.  Jobin)  performed  work  with  the  PI  and  Co-PI.  These  students 
were  entirely  supported  by  ONR  Award  #  N000141 512158.  Four  other  students  were  involved: 
Zachary  J.  Marinelli,  Aidan  L.  Gordon,  Brian  N.  Lang,  and  Yosyp  Schwab.  In  Fall  2015  Aidan 
L.  Gordon,  performed  modeling  and  data  analysis  work  for  this  Project  and  received  a  stipend 
for  it.  Supervised  by  the  PI,  Co-PI,  and  technical  personnel  at  JMU,  the  students  assembled  the 
insulated  compartment,  interfaced  the  experimental  set-up  with  the  electronic  data  collection 
system,  and  finally  collected  and  analyzed  the  experimental  data.  The  students  used  the  machine 
shops  at  JMU,  were  involved  in  the  intellectual  effort  of  modeling  and  understanding  the  data, 
and  presented  the  results  at  various  venues  [8-12],  The  student’s  efforts  were  rewarded  by  their 
inclusion  in  the  list  of  authors  of  the  publication  reporting  most  of  the  recent  findings  [4],  Five 
students  actively  participated  to  the  Workshop  on  “ Infrared  radiation,  Thermoelectricity  and 
Chaos ”  on  June  17,  2015  at  JMU. 

(•)  Academic  outreach :  In  Fall  2015,  the  PI  taught  the  course  on  Thermodynamics  and 
Statistical  Mechanics  PHYS  380  at  JMU  stemming  out  from  the  observation  that  the  conjugated 
variable-based  approach  can  (a)  broaden  the  concept  of  heat,  and  (b)  include  the  description  of 
several  phenomena  such  as  equilibrium,  stability,  the  laws  of  Thermodynamics,  phase 
equilibrium,  specific  heat,  compressibility,  boson  and  fennion  statistics,  Maxwell-Boltzman 
distribution  for  the  momentum  of  many  non-correlated  particles,  momentum  for  correlated 
particles,  the  law  of  ideal  gases,  properties  of  solids  at  low  temperatures,  magnetic  properties  of 
materials,  continuity  equation  in  fluids  and  Fick’s  laws.  The  PI  became  convinced  that  the 
conjugated  variable-based  approach  can  also  serve  as  the  framework  in  which  to  model  the 
research  in  this  Project  on  IR  power  generation.  The  conjugated  variable-based  approach 
generated  a  new  and  original  way  to  design  and  deliver  the  lectures  in  PHYS  380.  All  the  newly 
developed  class  material  is  now  in  hand-written  notes.  However  the  PI  thinks  that  it  is  worth  the 
effort  to  draft  a  textbook  based  on  these  notes  and  will  pursue  this  goal  in  the  future  months. 
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(•)  Broader  outreach;.  Two  main  objectives  were  achieved: 

— »  The  PI  illustrated  some  of  the  findings  of  the  research  project-ONR  Aw'ard  # 
N000141512158  in  an  event  organized  by  the  Mount  Jackson  Rotary  Club,  Mount  Jackson  (VA), 
on  August  4,  2015  [5].  Similar  events  might  take  place  in  the  near  future.  The  Pi’s  presentation 
was  well  received  and  effective  in  disseminating  the  awareness  of  the  U.S.  Navy’s  interest  in 
alternative  sustainable  energy  programs. 

— >  The  PI  presented  the  finding  in  the  project  at  two  international  venues  [6,  7]  focusing  on 
the  general  aspects  of  IR  power  generation  [6]  and  on  the  potentials  of  the  nano-PG  devices  [7]. 
The  Pi’s  presentations  generated  discussions  with  other  researchers  (e.g.  A.  V.  Zenkevich,  the 
organizer  of  the  event  in  Moscow)  and  created  connections  with  Angela  Belcher  -  MIT  in 
Boston,  Leora  Peltz  -  Boeing  in  Huntington  Beach  (CA),  and  Scott  Strum  -  Medical  School  at 
Loma  Linda  University  in  Loma  Linda,  CA. 
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THE  WORKSHOP  “ Infrared  radiation,  Thermoelectricity  and  Chaos”. 

The  Workshop,  which  took  place  on  June  17,  2015  at  JMU,  was  one  of  the  most  important 
activities  of  the  research  project-ONR  awards  #  N000141410378  and  N000141512158.  The 
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event  was  organized  by  the  PI  and  Co-PI  together  with  Mr.  Benjamin  T.  Delp,  Director  of  the 
JMU  Research  Development  and  Promotion  Office.  The  Workshop  featured  two  international 
speakers.  Prof.  Giulio  Casati  of  Italy  spoke  about  Conservation  laws,  symmetry’  breaking  and 
control  of  the  heat  current.  Prof.  Francisco  Javier  Gonzalez  Contreras  of  Mexico  described 
Seebeck  nanoantennas  for  solar  energy  harvesting.  Three  speakers  represented  the  U.S.A. 
academic  world.  Prof.  Gabriel  A.  Rincon-Mora  of  the  Georgia  Institute  of  Technology, 
presented  Microwatt  CMOS  harvesters.  Prof.  Rajeevan  Amirtharajah  of  the  University  of 
California  at  Davis,  discussed  Powering  systems  from  ambient  energy  sources.  Finally,  Prof  S. 
Joseph  Poon  of  the  University  of  Virginia,  presented  Practical  realization  of  ZT>1  from  the 
materials  perspective.  The  business  and  entrepreneurial  world  was  represented  by  Dr.  Ronald  J. 
Parise  of  Parise  Research  Technologies,  who  introduced  The  Nighttime  solar  cell' :  infrared 
energy  to  deep  space.  Dr.  Joseph  R.  Blandino  of  the  Virginia  Military  Institute,  and  Dr.  David  J. 
Lawrence  of  JMU,  presented  a  poster  on  Transient  response  of  a  thermoelectric  generator 
subjected  to  spatially  non  -uniform  heating:  implications  for  heat  and  IR  sensing  applications.  In 
her  presentation,  the  Dr.  G.  Scarel  (PI)  explained  the  connection  among  the  words  in  the 
Workshop’s  title  (IR  radiation,  thermoelectricity  and  chaos)  justified  the  selection  of  the 
speakers,  and  illustrated  the  current  knowledge  on  IRPG.  Mr.  Lynn  J.  Petersen,  ONR’s  Program 
Officer  (Code:331)  who  funded  the  awards,  participated  to  the  event. 

Two  posters  were  presented  by  JMU  Physics  students  Graham  P.  Gearhart  and  Brian  N. 
Lang.  Graham  presented  “ Infrared  source  and  thermoelectric  device  trends  in  infrared  power 
generation ”,  while  Brian  discussed  “ The  interaction  between  current  and  infrared  radiation”. - 
The  JMU  Physics  students  Justin  M.  Kaczmar,  Zachary  J.  Marinelli  and  Tara  R.  Jobin  joined  the 
event  by  helping  in  the  logistics,  asking  questions  during  the  sessions  and  actively  networking. 

Dr.  Brian  C.  Utter  (Co-PI)  and  Mr.  Lynn  J,  Petersen  (ONR)  chaired  the  morning  and  the 
afternoon  sessions,  respectively.  Dr.  Yvonne  R.  Harris,  JMU’s  vice  Provost  for  Research  and 
Scholarship,  and  Mr.  Lynn  Petersen  opened  the  event  in  the  morning  and  addressed  the  final 
remarks  in  the  afternoon,  respectively. 

The  Workshop  was  funded  by  the  ONR  award  #  N0001 4141 0378.  Additional  funding,  used 
to  support  the  logistics,  pay  the  honorarium  for  Dr.  J.  R.  Blandino  and  partially  that  of  Prof.  F.  J. 
Gonzalez,  was  provided  by  the  JMU  4-VA  Consortium  and  The  JMU  Madison  Trust — Fostering 
Innovation  and  Strategic  Philanthropy.  The  international  speakers  (Prof.  G.  Casati  and  Prof.  F.  J. 
Gonzalez)  were  mainly  funded  by  ONR  Global-London  (UK). 

The  Workshop  revealed  that  IR  power  generation  explores  a  rather  unknown  topic  in  a 
unique  experimental  setting.  Specifically: 

(•)  F.  J.  Gonzalez:  excites  a  PG  device  exploiting  Joule  heating,  i.e.  an  entropic 
contribution,  derived  from  an  electric  current  activated  in  an  antenna  by  the  IR  radiation.  His 
research  does  not  explore  the  direct  effects  of  radiation,  i.e.  of  the  electric  contribution,  on 
voltage  output  as  we  do  in  IR  power  generation. 

(•)  R.  J.  Parise:  uses  IR  radiation  as  a  means  to  carry  heat  to  deep  space  which,  acting  as  a 
heat  sink,  facilitates  the  creation  of  a  temperature  difference  in  a  PG  device. 

(•)  D.  J.  Lawrence  and  J.  R.  Blandino:  use  high  power  fiber  optics  to  activate  a  PG  device. 
Their  results  are  above  the  power  range  of  IR  pow'er  generation. 

(•)  G.  Casati:  pointed  out  that  the  term  “soliton”  was  improperly  used  by  PI  in  her 
presentation,  thus  underlining  that  the  hyperbolic  instability  highlighted  in  the  briefing  on  IR 
power  generation  is  a  new  entity. 
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FINANCIAL  REPORT. 

Expenditures. 

(•)  Stipends: 

— »  Giovanna  Scarel  (PI,  Assistant  Professor,  JMU,  Department  of  Physics  and  Astronomy), 
Summer  2015,  $  3,609.00 

— >  Brian  C.  Utter  (Co-PI,  Associate  Professor,  JMU,  Department  of  Physics  and 
Astronomy),  Summer  2015,  $  3,915.00 

— >  Justin  M.  Kaczmar  (Undergraduate  student,  JMU,  Department  of  Physics  and 
Astronomy),  May  21  —  July  31,  2015,  $  8.50/hour,  S  4,000.00 

— >  Graham  P.  Gearhart  (Undergraduate  student,  JMU,  Department  of  Physics  and 
Astronomy),  May  21  -  July  3,  2015,  S  8.50/hour,  S  2,000.00 

— >  Tara  R.  Jobin  (Undergraduate  student,  JMU,  Department  of  Physics  and  Astronomy), 
May  21  -  July  3,  2015,  $  8.50/hour,  S  2,000.00 

— »  Aidan  L.  Gordon  (Undergraduate  student,  JMU,  Department  of  Mathematics  and 
Statistics),  Fall  2015,  $  8.50/hour,  $  1,000.00 
(•)  Taxes: 

FICA:  $  1,187.59 
(•)  Travel: 

— >  Car  rental  from  JMU  for  G.  Scarel’s  travel  to  NRL  on  August  17-19,  2015:  $  82.42. 
(Sample  preparation  with  Virginia  Wheeler  and  data  collection  with  Marc  Currie). 

— >  Lodging  and  food  for  G.  Scarel’s  travel  to  NRL  on  August  17-19,  2015:  $  655.82 
(•)  Miscellaneous: 

— >  Bruker  HeNe  laser  module  to  support  IR  source  in  Pi’s  laboratory;  November  2015:  $ 
1,500.00 

— »  PIKE  Technologies  Crystal  Polishing  and  Polishing  kit  Cerium  Oxide  compound  to  clean 
KBr  window  in  IR  source  compartment  in  Pi’s  laboratory;  December  2015;  $  78.42 

— »  CDW-Govemment  AME-BE350G  APC  Back  UPS  for  power  support  of  Pi’s  IR  source  in 
sample  compartment;  December  2015:  $  54.01 

(•)  Indirect  cost: 

— >  $  9,173.00 
Other  support. 

(•)  4VA  Collaborative  Research  Project  2014-2015:  “ Giving  to  JMU  international  visibility 
in  infrared  power  generation  and  its  link  to  non-linear  phenomena” .  PI:  G.  Scarel  (JMU).  Co- 
PIs:  Rosa  A.  Lukaszew  (College  of  William&Mary)  and  Brian  C.  Utter  (JMU).  (“ Infrared 
radiation,  Thermoelectricity  and  Chaos"  Workshop  support). 

(•)  JMU  -  Madison  Trust — Fostering  Innovation  and  Strategic  Philanthropy  -  Innovation 
Grant  -  2015,  November  14,  2014:  “ From  JMU  to  the  World:  embracing  the  need  of  new  energy 
sources".  (Student  support  and  miscellaneous  expenses  for  the  Workshop  on  "Infrared  radiation, 
Thermoelectricity  and  Chaos"  Workshop). 

(•)  JMU  -  Department  of  Physics  and  Astronomy.  (Student,  equipment  and  facility  support). 

ATTACHMENTS. 

1)  Presentations  (3) 

2)  Papers  (1) 
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ABSTRACT. 

This  project  studies  the  basic  processes  involved  in  the  interaction  of  infrared  radiation  with 
matter  and,  in  this  effort,  establishes  a  method  to  harvest  infrared  radiation  and  transform  it  into 
usable  energy.  We  have  identified  an  electrical  contribution  to  the  energy  transfer  from  the 
infrared  radiation  to  the  device  we  use  to  harvest  it.  This  device  behaves  as  a  capacitor,  and  we 
found  a  way  to  triplicate  the  voltage  generated  by  the  device  when  exposed  to  the  infrared 
radiation.  This  significant  increase  can  be  obtained  by  placing  two  devices  in  series  and 
covering  the  surface  of  the  device  facing  the  infrared  radiation  to  tape,  which  modifies  the 
capacitance  of  the  device.  We  plan  on  modeling  the  interaction  of  the  infrared  radiation  with  our 
device  using  the  conjugated  variable-based  approach.  We  also  plan  on  studying  the  current 
produced  by  the  device  to  increase  efficiency  in  power  generation.  We  are  also  looking  into 
other  alternative  device  designs  to  fabricate  nano-  and  macro-power  generation  device.  Seven 
undergraduate  students  were  involved  in  the  research  and  the  research  inspired  an  approach  to 
Thermodynamics  and  Statistical  Mechanics  related  to  the  conjugated  variable-based  approach. 
We  plan  on  outlining  this  view  in  a  textbook. 

Keywords:  infrared  radiation,  radiation  harvesting,  renewable  and  sustainable  energy  source 
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Correlation  Between  the  Sinusoidal  Instability  in  Radiation  Power  and 
the  Hyperbolic  Instability  in  Voltage  in  Infrared  Power  Generation 

Justin  Michael  Kaczmar1  Tara  Rose  Jobin1,  Brian  C.  Utter1  and  Giovanna 

Scarel1 

’Department  of  Physics  and  Astronomy,  James  Madison  University, 

Harrisonburg  VA  22807 


Abstract: 


Goal-  Looking  for  methods  to  harvest 


alternative  sources  of  energy,  specifically 
infrared,  and  transform  it  into  usable 
energy.  The  advantage  in  harvesting 
infrared  is  its  day  and  night  availability. 
Specifically  we  are  looking  at  the 
instability  in  the  power  source  in 
comparison  to  the  instability  in  the 
resulting  voltage  the  device  produces. 


22  55 


Experimental  setup: 


10  20  30  40 

Time  [h] 


Infrared  Radiation: 

Device  AIO  plate 


Observing  a  sinusoidal, 
instability  in  the  power 
source. 


similar  tc  times 


Equation:  Voltage 


AF(0  =  AF/2+AF05C  *sec//^j 

I 


Results: 


10  20  30 

Time  [h] 


Power 


Voltage 


*tc  is  the  critical  time 
of  both  of  the  systems 
according  to  their 
respective  functions 
where  they  are  at  a 
minimum 


0.972 


There  is  a  correlation  between  the 
parameters  between  the  two  functions, 
specifically  with  tc. 
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Abstract: 


The  nano-power  generator  fabricated 
with  thin  atomic  layer  deposited 

G.  Scarel.1  V.D  Wheeler,2  H.S  Mann,1  B.N  Lang,1  Z  J  Marinelli,1  and  B.C. 

Utter1 
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Experimental  setup: 


Can  infrared  (IR)  radiation  [1,2] 
generate  voltage  on  a  nano-power 
generator  (NPG)  with  active 
element  consisting  of  an  atomic- 
layer  deposited  (ALD)  film  [3,4]? 


Anodized  aluminum 
sample  holder 


ALD  Nb-doped  DO, 
thermoelectric  film 


Nanometric  power 
■  generator  (NPG) 
device 


IR/Hoat 


Result-1 : 


Result-2: 
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Result-3: 


68.4  nm  thick  Nb-doped  Ti02  on 
borosilicate  glass, 

resistivity  =  1.4  10-5  Q  m  (~  graphite) 


0  005 


-0.005 


Timers] 


1000 


H.S.  Mann  etal.,  J.  Vac.  Sci.  Technol.  A  33,  01A124 
(2015). _ 


Ti02:TiN  ~2  nnr-2  nm  on  glass, 
sheet  resistance  =  26700  Ohms/sq 
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These  preliminary  results 
suggest  that  power 
generation  with  ALD 
NPG  might  be 
significantly  improved 
by  engineering  the 
film  resistance. 


I.  The  IR  radiation  is  able  to  excite  the 
voltage  on  a  NPG. 

II.  Voltage  and  temperature  difference 
have  the  same  time  constant. 

III.  Correct  resistivity  to  achieve  larger 
voltage  rise  upon  starting  the 
illumination  with  IR  radiation  still 
needs  to  be  found. 
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Abstract 


The  interaction  between  infrared  radiation  and  a  power  generator  device  in  time  is  studied  as  a 
route  to  harvest  infrared,  and  possibly  other  electromagnetic  radiations.  Broadening  the  spec¬ 
trum  of  the  usable  electromagnetic  spectrum  would  greatly  contribute  to  the  renewable  and  sus¬ 
tainable  energy  sources  available  to  humankind.  In  particular,  low  frequency  and  low  power  radi¬ 
ation  is  important  for  applications  on  ships,  satellites,  cars,  personal  backpacks,  and,  more  generally, 
where  non-dangerous  energy  is  needed  at  all  hours  of  the  day,  independent  of  weather  conditions. 
In  this  work,  we  identify  an  electric  and  an  entropic  contribution  to  the  energy  transfer  from  low 
power  infrared  radiation  to  the  power  generator  device,  representing  electrical  and  thermal  con¬ 
tributions  to  the  power  generation.  The  electric  contribution  prevails,  and  is  important  because  it 
offers  multiple  ways  to  increase  the  voltage  produced.  For  example,  placing  black-colored  gaffer 
tape  on  the  illuminated  face  doubles  the  voltage  produced,  while  the  temperature  difference,  thus 
the  entropic  contribution,  is  not  sensitive  to  the  presence  of  the  tape.  We  recognize  the  electric 
contribution  through  the  fast  changes  it  imparts  to  the  voltage  output  of  the  power  generator  de¬ 
vice,  which  mirror  the  instabilities  in  time  of  the  infrared  radiation.  The  device  thus  acts  as  sensor 
of  the  infrared  radiation's  behavior  in  time.  On  the  other  hand,  we  distinguish  the  entropic  con¬ 
tribution  through  the  slow  changes  it  causes  to  the  voltage  output  of  the  power  generator  device, 
which  reflect  the  relative  delay  with  which  the  two  faces  of  the  device  respond  to  thermal  pertur¬ 
bations. 
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1.  Introduction 

A  power  generator  (PG)  device  can  be  used  to  harvest  electromagnetic  (EM)  and,  in  particular,  infrared  (IR) 
radiation.  The  interaction  between  the  radiation  and  the  device  is  a  complex  phenomenon  of  energy  transfer 
(AE).  The  rate  of  energy  transferred  from  the  EM  radiation  per  area  a  of  the  device  is  the  Poynting  vector 
|£|  =  |£x//|  =  P/a,  where  E  and  H  are  the  electric  and  magnetic  fields,  respectively,  and  P  is  power. 
Therefore,  because  of  E  and  //  ,  the  interaction  between  radiation  and  device  involves  the  charges  on  the  de¬ 
vice  surface.  Electromagnetic  radiation  with  large  frequency  v  interacts  through,  e.g.,  Compton  scattering  [I], 
X-ray  photoelectron  effect  [2],  photoelectric  effect  [3],  photovoltaic  effect  [4],  and  plasmon  generation  [5]. 
Electromagnetic  radiation  with  low  frequency  v ,  e.g.  in  the  IR  and  microwave  regions,  resonates  with  molecu¬ 
lar  rotation  and  oscillation  frequencies  [6]  or  generates  polaritons  [7]-[10].  When  the  photon  frequency  v  or 
energy  hv ,  where  h  is  Planck’s  constant,  do  not  match  with  the  frequency  or  the  energy  of  a  specific  pheno¬ 
menon  involving  charges,  the  energy  of  the  EM  radiation  contributes  to  temperature  T  changes.  In  photosynthe¬ 
sis  this  phenomenon  is  known  as  internal  conversion  [11]. 

We  name  the  energy  transferred  from  the  EM  radiation  to  a  PG  device  through  the  action  of  the  electric  E 
and  magnetic  //  fields  as  the  electric  contribution: 

AEel  -qAV ,  (1) 

where  q  is  the  charge  and  V  voltage.  We  name  the  energy  transferred  through  changes  in  temperature  T  at  en¬ 
tropy  X  as  the  entropic  contribution: 

AEen=ZA  T.  (2) 

The  energy  transferred  from  IR  and  microwave  radiation  is  usually  associated  with  the  entropic  contribution 
in  Equation  (2).  For  example,  sun  light  gives  the  sensation  of  temperature  increase,  and  therefore  of  warmth,  on 
human  skin.  The  microwave  radiation  in  microwave  ovens  is  used  to  increase  the  temperature,  i.e.,  cook  food 
and  heat-up  beverages.  Similarly,  through  laser  radiation  it  is  possible  to  increase  temperature,  even  with  na¬ 
noscale  control  [12]. 

The  effects  of  the  electric  contribution  AEel  -qAV  are  less  apparent  in  the  energy  transfer  from  low  fre¬ 
quency  and  low  power  EM  radiation.  In  the  current  literature,  the  existence  of  the  electric  contribution  is  ac¬ 
knowledged  [  1 3]-[  1 6],  but  the  interplay  between  the  electric  and  the  entropic  contributions  is  not  investigated. 
Specifically,  there  is  a  lack  of  knowledge  of  1)  the  possibility  of  decoupling  the  electric  from  the  entropic  con¬ 
tributions,  2)  the  factors  that  promote  the  electric  over  the  entropic  contribution,  or  vice-versa ,  3)  the  existence 
of  a  threshold  where  one  contribution  prevails  over  the  other,  and  4)  the  benefits  of  the  electric  over  the  entropic 
contributions,  or  vice-versa. 

In  this  work  we  aim  at  decoupling  AEel  and  AEen  in  a  PG  device  illuminated  by  low  power  IR  radiation. 
The  device  is  expected  to  respond  to  the  entropic  contribution  by  exploiting  the  Seebeck  effect  [17]-[20],  i.e. 
producing  a  voltage  difference  AV  directly  proportional  to  the  temperature  difference  AT  applied  to  the  two 
faces  of  the  PG  device,  so  that  AV  = -SAT .  Here,  S  is  the  Seebeck  coefficient.  On  the  other  hand,  we  expect 
the  PG  device  to  also  respond  to  the  electric  contribution  through  its  capacitor-type  of  structure  consisting  of  a 
sequence  of  conducting  and  insulating  layers,  as  illustrated  in  Figure  1.  For  the  device  used  in  this  work,  the 
sequence  is,  starting  from  the  face  illuminated  by  the  IR  radiation,  a  copper  (Cu)  plate,  a  layer  of  pillars  made  of 
adoped  Bi2Te3-based  alloy,  another  Cu  plate,  and,  finally,  an  alumina  (AlO)  plate.  On  the  Cu  plates  there  are 
electrons  whose  surface  density  a=q/a  is  sensitive  to  the  E  and  II  fields  of  the  IR  radiation,  thus 
enabling  changes  in  the  electric  contribution  AEel  =  qAV . 

In  our  experiment,  the  voltage  difference  AE(/)  ,  generated  by  the  PG  device  through  the  electric  and  the 
entropic  contributions,  and  the  temperature  difference  AT(/)  ,  related  to  the  entropic  contribution,  are  observed 
as  a  function  of  time  /.  The  measurements  capture  the  first  minutes  after  starting  the  illumination,  and  in  the  30 
hours  thereafter.  We  hypothesize  that  changes  in  |5(/)j  x  a  =  |£(/)x  //(/)|x  a  =  P(/)  slowly  vary  the  ampli¬ 
tude  of  the  surface  electron  density  cr(r,/).  To  prove  this  hypothesis,  we  study  the  power  P(/)  of  the  IR 
radiation  using  a  power-meter  sensor  and  compare  its  behavior  with  that  of  AV  (/)  and  AT(t)  . 

Summarizing,  we  consider  the  total  energy  transfer  in  time  from  the  IR  radiation  to  a  PG  device  as  the  sum  of 
the  electric  and  the  entropic  contributions  such  that: 


(3) 
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(a)  <b> 


Figure  1.  Schematics  of  the  away  (a)  and  toward  (b)  architectures  of  the  PG  device  In  the  away  architecture  (a)  the  face  of 
the  PG  device  exposed  to  the  IR  radiation  is  free  from  contact  with  the  sample  holder.  In  the  toward  architecture  (b),  the  il¬ 
luminated  face  is  in  contact  with  the  sample  holder.  The  PG  device  is  a  stack  of  conducting  (Cu  plates),  non-conducting 
(AlO  plate),  and  semiconducting  (set  of  pillars  made  of  a  doped  Bi2Te3-based  alloy)  layers. 


Consequently,  we  assume  the  voltage  difference  AV  (V) 
tion  of  two  summands: 


produced  by  the  PG  device  in  time  to  be  the  addi- 


AF(/) 


*(') 


ZAT(t)' 

~wr) 


(4) 


The  first  summand  relates  to  the  electric  and  the  second  to  the  entropic  contribution.  The  term  S  =  —Y,jq(t) 
can  be  associated  with  the  Seebeck  coefficient 

We  will  show  that  with  the  low  power  irradiation  employed  in  our  measurements,  the  electric  contribution  can 
be  decoupled  from  the  entropic  contribution,  and  largely  dominates.  Decoupling  the  two  contributions  is  impor¬ 
tant  for  IR  energy  harvesting,  because  the  electric  contribution  offers  a  variety  of  ways  :o  increase  the  voltage 
produced  by  the  PG  device,  e.g.  by  placing  black-colored  gaffer  tape  on  the  illuminated  face  of  the  device,  as  we 
will  show  in  Appendix-1.  The  entropic  contribution,  nstead,  is  limited  by  the  temperature  difference  AT(/) 
established  between  the  two  faces  of  the  PG  device. 


2.  Experimental  Set-Up 

For  this  experiment,  continuous  broadband  IR  radiation  in  the  middle  IR  (MIR)  region  (i.e.  frequency  between 
350 -7500  cm  1 ,  or  wavelength  between  20 -2.2 pm)  was  produced  by  a  globar  (Q301)  source.  The  power 
of  the  IR  radiation  was  monitored  versus  time  udng  a  power-meter  sensor  Coherent  Power  Max  RS  PS  19, 
sensitive  to  the  300  - 1 1 000  nm  wavelength  range,  and  to  the  100  pW  to  1  W  power  range. 

The  voltage  difference  A V (t) ,  generated  by  the  electric  and  the  entropic  contributions  to  AEtot  according 
to  Equation  (4),  was  produced  using  a  PG  device  071 1 1-9L31-04B  by  Custom  Thermoelectric  Inc.  The  device 
consists  of  a  sequence  of  layers:  1)  a  Cu  plate  on  the  face  exposed  to  the  IR  radiation,  2)  a  layer  of  pillars  made 
of  a  doped  Bi2Te3-based  alloy,  3)  another  Cu  plate,  and  4)  an  AlO  plate.  The  Cu  plate  no:  illuminated  by  the  IR 
radiation  is  non-continuous,  as  highlighted  through  the  white  hole  in  the  left  side  of  Figure  1(a)  and  Figure  1(b). 
In  the  away  architecture,  illustrated  in  Figure  1(a),  we  established  the  continuity  by  placing  the  sample  holders, 
made  of  anodized  aluminum,  in  contact  with  the  non-continuous  Cu  plate.  Thus,  the  Cu  plate  together  with  the 
sample  holder  behaves  as  the  electrode  of  a  capacitor,  ^he  illuminated  Cu  plate,  instead,  was  free  of  contact  with 
the  sample  holder.  In  the  toward  architecture,  pictured  in  Figure  1(b),  we  left  non-continuous  the  Cu  plate  op¬ 
posite  to  the  IR  radiation,  while  the  illuminated  Cu  plate  was  kept  in  contact  with  the  sample  holder. 

The  temperatures  TiR  (/)  and  Tno_m(t)  of  the  illuminated  and  non-ilium  ’mated  faces,  respectively,  of  the 
PG  device  were  measured  using  OMEGA  type  E  Ni-Cr/Cu-Ni  thermocouple  probes.  The  temperature  difference 
A  T(t)  was  obtained  as  AT (/)  =  Tm  (/)  -  Tno  IR  (/)  .  T ae  trends  of  AV (/) ,  Tm  (/)  and  Tno_m  (7)  were  meas¬ 
ured  using  Keithley  2000  multi-meters.  The  data  were  collected  using  LabView  2012  and  a  National  Instru¬ 
ments  PXI-1042q  communications  chassis. 

During  the  measurements,  the  PG  device  and  the  power-meter  sensor  were  positioned  vertically  and  at  an  an¬ 
gle  of  incidence  0Q  =  45  with  respect  to  the  IR  radiation.  The  instrumentation  was  placed  in  a  closed  sample 
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compartment  purged  with  N2  to  prevent  disturbances  for  the  whole  duration  of  the  measurements  versus  time  of 
H‘)'  A V{t),  Tir  (/)  and  Tno_m  (/“)  [21].  The  experimental  parameters  are  summarized  in  Table  1. 

3.  Results  and  Discussion 

a)  Behavior  in  time  of  P(t) 

In  the  100  seconds  immediately  after  starting  the  illumination  of  the  power-meter  sensor,  P(t)  ,  displayed  in 
Figure  2(a),  rises  exponentially  as  follows: 


Tabic  1.  Summary  of  the  experimental  parameters  in  the  main  text  and  in  the  appendices 


Main  text 

IR  frequency  (pm) 

Infrared  angle  of  incidence  (  ) 

PG  device  final  layer 

20  -  2.2 

45 

Cu 

Appendices 

IR  frequency  (pm) 

Infrared  angle  of  incidence  (  ) 

PG  device  final  layer 

20-2.2 

30 

Alumina 

Figure  2.  (a)  Exponential  rise,  as  in  Equation  (5),  of  the  power  P(t )  versus  time  of  the  IR  radiation  emitted  by  the  globar 
source  in  the  100  seconds  immediately  after  starting  the  illumination  of  the  power-meter  sensor;  (b)  Graph  of 


d  / 


in  the  same  time  interval  of  (a)  reporting  the  slope  and  amplitude  Ap ;  (c)  The  power  7>(/)  in  the  50 


hours  after  starting  the  illumination  of  the  power-meter  sensor;  (d)  Same  as  (c),  with  the  vertical  scale  expanded  to  highlight 
the  sinusoidal  instability  region  fitted  with  Equation  (6).  The  zero  of  the  time-scale  coincides  with  the  start  of  the  illumina¬ 
tion  with  the  IR  radiation.  The  parameters  tcF  and  Up  are  labeled;  (e)  Graph  of  \P(t),^—^—\  in  the  50  hours  after 
starting  the  illumination  of  the  power-meter  sensor 
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P(t)=P0e^+Pf, 


(5) 


where  P0  is  the  offset,  Pf  the  final  value,  and  rp  the  time  constant.  Typical  values  of  these  parameters  are 


reported  in  Table  2.  The  slope  of  the 


p0b 


dp(0 

di 


graph  in  Figure  2(b)  is  negative,  indicating  the  evolution 


of  P(t )  toward  a  stable  fixed  point  [22].  The  rate  of  increase  of  ,  i.e.  the  absolute  value  of  the  slope  of 

(  ,  dp(iY 

n*)-— 


the 


d  t 


graph,  is  pp=  2.04  s  1 .  The  amplitude  Ap  of  P(t ),  i.e.  the  magnitude  of  the  interval 


along  the  horizontal  scale  of  the 


Hi) 


0,^ 


d  / 


graph  in  Figure  2(b),  is  21.5  mW.  In  the  50  hours  after  start¬ 


ing  the  illumination  of  the  power-meter  sensor,  the  power  P(/) ,  shown  Figure  2(c),  reaches  a  plateau.  Howev¬ 
er,  in  multiple  data  sets,  we  always  observe  that  /*(/)  undergoes  small  sinusoidal  instabilities  shown  in  Figure 
2(d)  in  which  the  vertical  scale  has  been  expanded.  The  instability,  due  to  small  periodic  fluctuations  in  the 
closed  sample  compartment,  can  be  fitted  with: 


p{‘)=  Poff  +  PoscS  in 


Hp  J 


(6) 


where  Po//  and  Posc  are  the  offset  value  and  half  of  the  separation  between  PMAX  and  Pmm ,  respectively. 

The  critical  time  tcP  is  the  point  in  time  in  which  p(t)  reaches  Pmn  +  Posc ,  while  Hp  is  the  amount  of  time 

necessary  to  move  P(t)  from  Pmn  +  Posc  to  EMAX  (practically,  1/4  of  the  period  of  the  sinusoidal  function). 
The  typical  values  of  Poff  ,  POK  ,  tcP  ,  and  Hp  are  reported  in  Table  2,  labelled  in  Figure  2(d),  and  were  ob¬ 
tained  by  placing  the  zero  of  the  time-scale  at  the  start  of  the  illumination.  The  instability  in  P(t)  is  small,  as 

(  An(  *W 


inferred  from  Poff  / Posc  »  103 .  The 


graph,  shown  in  Figure  2(e),  highlights  a  periodic  behavior 

V  al  J 

with  frequency  vp  =  \/AH p  =  2x10  5  Hz  !  We  observ'e  the  sinusoidal  instability  of  P{t) ,  which  modulates  the 
amplitudes  of  the  electric  E  and  magnetic  H  fields  of  the  IR  radiation,  to  persist  beyond  the  50  h  time  in¬ 
terval  in  Figure  2(d). 

b)  Action  of  the  IR  radiation  on  a  (/*,/) 

We  observed  that  the  power  P{t )  of  the  IR  radiation  rises  exponentially  obeying  Equation  (5)  at  the  start  of 
the  illumination,  and  exhibits  a  sinusoidal  instability  in  the  50  hours  thereafter.  For  the  entire  time  span,  we  hy¬ 
pothesize  that  the  IR  radiation  transfers  energy,  through  electric  contribution  AEel,  to  the  surface  density 


Table  2.  (Top  rows)  Fitting  parameters  P0 ,  Pf ,  and  rP  of  the  IR  power  =  P0e  Tp  +  Pf  in  Equation  (5)  in  the  100 

seconds  immediately  following  the  start  of  the  illumination  of  the  power-meter  sensor  with  IR  radiation  from  the  globar 
source.  The  rate  of  increase  of  the  power  (  pP  ),  and  the  amplitude  Ap  in  this  time  interval,  derived  from  the 


£( o 

dt 


graph  in  Figure  2(b),  are  also  reported.  (Bottom  rows)  Typical  values  of  the  Poff ,  Posc ,  tcP  ,  and  HP  of 


the  sinusoidal  instability  P  (t )  =  PQjf  +  Posc  sin  j  bi  Equation  (6)  of  the  IR  power  P  (f )  in  the  50  hours  following  the 

start  of  the  illumination  of  the  power-meter  sensor  with  IR  radiation  from  the  globar  source.  The  parameters  tcP  and  Hp 
are  labelled  in  Figure  2(d).  The  unit  for  time  in  the  bottom  rows  of  this  table  is  the  hour  ([h]). 


Time  interval 

/“.[mW] 

Pf  [mW] 

4[mW] 

0-  100  s 

-21.9 

21.9 

0.5 

2.04 

21.5 

Time  interval 

pA  mW] 

/“..[raW] 

HF[h] 

‘A  h] 

- 

0  -  50  h 

21.7±0.7 

0.05  ±0.03 

4.0±0  5 

16±2 

- 
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cr(r,/)  of  the  electrons  on  the  illuminated  Cu  plate  of  the  PG  device,  and  contributes  to  producing  A V{t) 
through  the  E  and  H  fields  (electric  contribution)  and  Ar(/)  (entropic  contribution).To  prove  that  a  link 
exists  between  P(/)  and  AK(/),  possibly  also  between  P{t )  and  AT^),  we  sketch  the  behavior  of 
cr(r,/)  and  relate  it  to  the  observed  /*(/),  AE(/),and  AT(/). 

To  sketch  o-(r,/) ,  we  hypothesize  that,  while  hitting  the  surface  of  the  Cu  plate,  the  IR  radiation  modulates 
the  electric  field  E  through  the  sinusoidal  instability  of  the  IR  power  P(t)  =  a|S(/)|  =  a|£'(/)x  // (/)| .  In 

turn,  E  and  its  modulation  act  on  the  electrons  of  the  Cu  plate  with  force  f  =  eE ,  where  e  is  the  electron’s 
charge.  As  in  the  photoelectric  effect  [3],  /  displaces  the  electrons  away  from  the  location  in  which  the  IR 
radiation  impinges  on  the  Cu  plate,  locally  decreasing  their  surface  density'  such  that  a  (r,/)  oc  1  j />(/) .  However, 
unlike  in  the  photoelectric  effect,  /  does  not  kick  the  electrons  out  of  the  Cu  plate.  In  this  process,  cr^,/) 
varies  in  time  t  as  well  as  in  space  r ,  i.e.  the  2-dimensional  (2D)  surface  of  the  Cu  plate.  To  allow  us  versatility 
in  choosing  reference  system,  orientation  and  phase,  we  represent  the  2D  space  variable  r  as  the  complex  va¬ 
riable  z  =  rx  +ir  ,  where  /  is  the  imaginary  unit.  This  choice  resembles  that  adopted  to  describe  light  polariza¬ 


tion  through  Jones  matrices  [23]-[25].  Thus,  z  = 


V 

- 1 

X 

n> 

ra 

+irb 

_ry_ 

i  ^ 

a 

i _ 

Jc 

.  All  possible  rotations  of  the  refer¬ 


ence  system,  phases,  and  positions  in  the  2D  plane  can  be  obtained  by  selecting  magnitude  and  sign  of  ra ,  rb  , 
rc ,  and  rd  . 

With  this  choice  of  z ,  upon  starting  the  illumination,  we  picture  a(z,/)  to  exponentially  decrease  accord¬ 


ing  to  a(zj)  =  RE 


rx+iry)k+<p 


,  where  <r0  and  crf  are  the  initial  and  final  surface  electron  den¬ 


sities,  ra  the  time  constant,  A  a  vector  with  units  of  inverse  length,  and  (p  an  arbitrary  phase.  We  note  that 
the  exponential  behavior  is  modulated  by  the  oscillatory  function  Qrk . 

In  the  subsequent  30  hours,  from  Equation  (6)  we  expect  cr(z,t)  to  undergo  a  slow  variation  in  time  such 

y 

,  where  the  sine  function  has  the  frequency  vp  =  \/4ffp  derived  in  Sec- 


that  cr(z,/)  x  \J P{t}  x  \  sin 


[zLi 
v  H, 


tion  3(a).  With  the  choice  of  z  discussed  above,  and  utilizing  the  laws  of  trigonometric  functions  for  complex 
variables,  we  obtain: 


oc  RE 


vJ-v 


1  1 

=  RE 

(  \ 

1 

1  sech(rv,/) 

h  J 

(sin(z,/)J 

sin(rx  +  iry,y 

sin(rx,/)coshty,/)  sinfc,/)  ^ 

r  \ 

vj  -  V  J„y 

x  x  crx 

(7) 


Here,  vx  and  vy  are  the  instability’s  propagation  velocities  along  the  x  and  y  directions;  Lx  and  Ly  are 
the  lengths  of  the  Cu  plate  along  x  and  y ;  finally,  tax  and  /  are  the  critical  times  of  the  surface  electron 
density’s  instability  along  *  and  y.  Considering  croft  and  cr0SC ,  the  equilibrium  electron  density  and  its  devia¬ 


tion  from  equilibrium,  respectively,  we  obtain  cr(z,t)  =  <J  +cr 


sech 


sin 


.  The  2D  space  variable 


Z  ,  therefore  unveils  a  hyperbolic  instability  in  cr^,/)  modulated  by  a  sine  function. 

While  requires  a  spatiotemporal  set  of  variables,  the  functions  AF(/),  T/R(t ),  Tno  IR  (/) ,  and 

A  T(t)  are  only  time-dependent.  To  decouple  Z  from  the  effects  of  cr(z,/)  on  AF(/),  T//?(/),  Tm_IR  (/) , 
and  AT  (/) ,  and  allow  time  t  to  be  the  only  effective  variable,  we  integrate  cr(z,/)  over  the  surface  area  a  of 

the  Cu  plate  as  J  cr(z,t)dz  .  Because  of  the  capacitor-ty  pe  structure  of  the  PG  device,  with  overall  capacitance 
a\  <r(z,/)dz 

C,  we  expect  A l/(/)=— — — - .  This  integration  causes  the  loss  of  correlation  between  the  phase  of 
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P(t)  and  AF(/).  In  the  30  hours  after  starting  the  illumination,  since  l/sin(rc,/)  >  1  where  |sin(r^,/)|  *  0 , 
from  Equation  (7)  we  expect  the  sech(rv,t)  function  to  determine  the  behavior  in  time  of  AF(/)  .  We  also 
envision  the  capacitor- type  structure  of  the  PG  device  to  affect  7^(/)  and  Tno_IR{t ),  while  leaving  AT(/) 
constant. 

Summarizing,  we  expect  A K(/)  to  obey  an  exponential  behavior  in  the  first  minutes  after  starting  the  illu¬ 
mination,  and  to  exhibit  a  hyperbolic-secant-type  of  instability  in  the  30  hours  thereafter,  with  no  phase  rela¬ 
tionship  with  P(/)  ,  no  periodic  behavior,  and  with  AT(/)  constant.  We  are  currently  exploring  this  hypothe¬ 
sis  further  numerically,  which  will  be  the  focus  of  future  work. 

c)  Behavior  of  AV (/)  and  AT(/)  in  the  400  seconds  after  starting  the  illumination 

The  voltage  difference  A K(/)  and  temperature  difference  AT(/)  in  the  400  seconds  after  starting  the  il¬ 
lumination  of  the  PG  device  with  1R  radiation  are  pictured  in  Figure  3  and  Figure  4  for  the  away  and  the  to¬ 
ward  architectures.  We  find  that  the  voltage  difference  AF(/)  in  Figure  3(a)  and  Figure  3(d)  fits  a  sum  of 
exponential  functions: 


*r(0-I 


AK0,e  r'  +  AF, 


(8) 


where  the  AK0/  and  AX'fi  terms  are  the  initial  and  final  voltage  differences,  and  the  zi  terms  are  the  time 
constants.  These  parameters  and  the  number  of  summands  N  are  summarized  in  Table  3.  In  both  the  away  and 
the  toward  architectures,  the  electric  contribution  is  related  to  the  summand  with  N  =  1  and  rx  =  3.8  s  .  The 


Figure  3.  Panels  (a),  (b),  and  (c)  correspond  to  the  away  architecture  and  refer  to  the  400  seconds  immediately  tollowing  the 
start  the  illumination  of  the  PG  device  with  1R  radiation,  (a)  Voltage  difference  AT(r)  with  fitting  curves  obeying  Equa¬ 
tion  (8)  highlighting  the  summands  related  to  the  electric  (el)  and  the  entropic  (en)  contributions;  (b)  Graph  of 
dA  v(ty\ 


AF(/),- 


dt  ) 


obtained  from  the  fitting  parameters  in  Table  3,  reporting  the  slope  and  amplitude  A;  (c)  Dimensionless 


voltage  AK(T)  =  AV(t)jAVfl  as  in  Equation  (9).  Panels  (d),  (e),  and  (f)  report  the  voltage  difference  A V(t)>  the 

AV graph  with  slope  and  amplitude  A,  and  the  dimensionless  voltage  AF(T)  =  AV{t)jAVfl ,  respectively, 

for  the  toward  architecture  in  the  400  seconds  immediately  following  the  start  the  illumination  of  the  PG  device  with  IR  rad¬ 
iation.  Panel  (d)  highlights  the  two  summands  related  to  the  electric  (el-1  and  el-2)  contributions,  and  the  summand  related  to 
the  entropic  (cn)  contribution. 
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entropic  contribution  is  related  to  the  summand  with  larger  time  constant,  r2  =  55  s  ,  and  is  labeled  with  N  —  2. 
This  association  is  justified  because  r2  =  55  s  is  of  the  same  order  of  magnitude  of  rn  =  30  s  for  the  temper¬ 
ature  difference  Af(/)  in  both  architectures,  as  can  be  seen  in  Table  3,  and  further  in  Appendix- 1.  The  sum¬ 
mand  with  N  =  3,  detected  only  in  the  toward  architecture,  is  related  to  an  electric  contribution  because  the 
“decay”  in  A  K(/)  ,  shown  in  Figure  3(d),  is  absent  in  the  corresponding  A  f(/)  in  Figure  4(c).  We  claim  that 
the  summand  with  N  =  3  causes  the  “decay”  in  A  V  {t)  shown  in  Figure  3(d)  because  of  the  non-continuity 
of  the  Cu  plate  opposite  to  the  1R  radiation  illustrated  in  Figure  1(b).  The  large  value  of  AF03  =  1.1  mV  signals 
that  the  toward  architecture  tries  to  avoid  the  “decay”  by  pumping  up  the  voltage  production 

c  d^  j 

3(b)  and  Figure  3(e),  are  p^v  -  0.25  s  1  in  the  away  architecture,  and  p^v  -  0.29  s  1  in  the  toward  architec¬ 
ture.  The  amplitudes  A ,  also  derived  from  Figure  3(b)  and  Figure  3(e),  are  0.95  mV  and  0.66  mV  in  the  away 
and  the  toward  architectures,  respectively. 


The  rates  of  increase  of  AV  (/)  ( p^v )  in  the  initial  linear  regime  of  the 


graph  in  Figure 


Tabic  3.  Fitting  parameters 


r,,  AF0|. ,  and  AVfi  of  the  voltage  difference  AT  (/)  in: 


^(0=1 


Afe  r'  +AV„ 


Eq¬ 


uation  (8)  (top  rows)  and  rTI ,  A Ton  and  A Tfi  of  the  temperature  difference  A T(t) 


A^(/)  =  Z 


A7;,e 


rn  +A7> 


Equation  (10)  (bottom  rows)  in  the  away  and  toward  architectures  in  the  400  seconds  immediately  following  the  start  of  the 
illumination  of  the  PG  device  with  1R  radiation.  The  corresponding  experimental  data  are  shown  in  Figure  3(a)  and  Figure 
3(d),  and  Figure  4(a)  and  Figure  4(c).  The  indexes  N  and  M  indicate  the  number  of  summands  in  Equations  (8)  and  (10), 
respectively.  The  relationship  of  the  summands  with  either  the  electric  or  the  entropic  contribution  is  highlighted. 


Away 

AV ;  N  =  2 

Figure  3(a)  jV  =  1  .  electric 

Figure  3(a);  N  —  2  ,  entropic 

* 

r,  -  3.8  s 

r2  =  55.0  s 

- 

AF0I  =  -0.75  mV 

A  Vn=  -0.20  mV 

- 

AV^  =0.78  mV 

AVf2  =0.30  mV 

- 

Toward 

AV  ;  N-3 

Figure  3(d):  N  -  \ ,  electric 

Figure  3(d)  N  =  2  ,  entropic 

Figure  3(d)  V  =  3  electric 

rt  =3  8s 

r2  =  55.0  s 

t3  =  48.0  s 
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- 

- 
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- 

- 
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- 

- 
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Figure  4.  Panels  (a)  and  (b)  correspond  to  the  away  architecture  and  refer  to  the  400  seconds  immediately  following  the  start 
of  the  illumination  of  the  PG  device  with  IR  radiation,  (a)  Temperature  difference  A7’(/)  with  exponential  behavior  as  in 

dArtOl 


Equation  (10);  (b)  Graph  of  AT(/) 


dt 


obtained,  from  fitting  parameters  in  Table  3,  reporting  the  slope  and  ampli¬ 


tude  Ar .  Panels  (c)  and  (d)  report  the  temperature  difference  A T(t)  and  the  |  A7'(/)?— • 


graph  with  slope  and 


amplitude  Ar  for  the  toward  architecture  in  the  400  seconds  immediately  following  the  start  the  illumination  of  the  PG  de¬ 
vice  with  1R  radiation.  Panel  (c)  highlights  the  two  summands  related  to  the  entropic  (en-1  and  en-2)  contribution. 


As  with  the  power  P(/) ,  the  voltage  difference  AF(/)  produced  by  the  PG  device  also  rises  exponentially 
in  the  first  few  minutes  after  starting  the  illumination.  However,  p^v-  0.25  and  0.29  s  r  are  one  order  of 
magnitude  smaller  than  pp  -  2.04  s~l  for  the  power  P(t)  shown  in  Figure  2(b).  We  ascribe  this  discrepancy 
to  the  activation  process  occurring  in  the  PG  device  to  start  the  AV  (/)  production. 

We  capture  the  relative  behavior  of  the  electric  and  the  entropic  contributions  to  the  generation  of  A F(/) 

using  the  dimensionless  voltage  AV {!)  =  AV (t) j AV f2  [22],  pictured  in  Figure  3(c)  and  Figure  3(0,  and  de¬ 
rived  from  Equation  (8)  and  its  parameters  in  Table  3  as  follows: 


AF(T)  =  t 
<=1 


vAP7; 


AV, 


/2 


(9) 


Here,  the  term  T  =  //r2  is  the  dimensionless  time  [22],  while  r0i  =  AF0,/AK/2  and  rfi=  AV rJ AV f2  are 
dimensionless  voltage  parameters.  We  performed  the  normalization  with  respect  to  the  slower  summand,  i.e.  the 

one  related  to  the  entropic  contribution  with  N  =  2  .  Such  a  choice  leads  us  to  observe  an  evolution  of  AF(T) 
in  Figure  3(c)  and  Figure  3(0  on  a  similar  T-scale  in  both  the  away  and  the  toward  architectures,  signifying 
that  the  dynamics  is  not  affected  by  the  architecture.  On  the  other  hand,  from  the  plateau  value  at  AF(T) «  4 

in  the  away  architecture,  we  infer  that  the  summands  related  to  the  electric  and  entropic  contributions  add  one  to 
the  other,  with  the  former  prevailing  over  the  latter  because  of  rn  -  AVfXjAVf2  =  2.6  >  1 .  On  the  contrary,  from 
the  plateau  value  at  AF(T)  <  1  in  the  toward  architecture,  we  conclude  that  the  two  summands  related  to  the 

electric  contributions  with  A-l  and  3  almost  annihilate  one  another,  while  the  summand  related  to  the  en¬ 
tropic  contribution  survives  in  the  long  term. 

The  temperature  difference  A7”(/),  shown  in  Figure  4(a)  and  Figure  4(c),  surprisingly  exhibits  negative 
values,  indicating  that  Tm  (/)  <  Tno_m  (/) .  This  behavior  is  typical  when  the  PG  device  faces  the  IR  radiation 
with  the  Cu  plate.  Nevertheless,  we  find  that  the  temperature  difference  AT  (/)  rises  exponentially  as  follows: 
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Ar(0-S 


A  T0lc  r-  +  A7^ 


(10) 


are  the  initial  and  final  temperature  differences,  respectively,  and  the  rTi  terms  are  the 

Tm  (/)  rises  faster  than  Tno_m  (/)  .  The  pa- 


vvhere  A  T0j  and  A Tf 

time  constants.  These  time  constants  reflect  the  fact  that  generally 
rameters  and  the  number  of  summands  A/,  all  related  to  the  entropic  contribution,  are  summarized  in  Table  3 
We  observe  that  I  <  M  <  2 ,  and  that  summands  similar  to  that  with  Af  =  2  (rarely  detected)  in  the  toward  ar¬ 
chitecture  in  Figure  4(c)  exhibit  an  extremely  large  time  constant  of  rT2  =  2000  s .  We  also  note  that  no  “decay” 
appears  for  AT(/)  in  the  toward  architecture  in  Figure  4(c),  as  opposed  to  the  corresponding  voltage  differ¬ 
ence  AK(/)  in  Figure  3(d). 

dA 7^)) 

d'  J 


The  slopes  of  the  initial  linear  regimes  in  the 


A?  (/)>“ 


graphs,  displayed  in  Figure  4(b)  and 


Figure  4(d),  are  negative,  as  those  of  /*(/)  in  Figure  2(b),  and  A K(/)  in  Figure  3(b)  and  Figure  3(e).  The 
rates  of  increase  of  AT(/)  are  pAr=  0.033  s1  and  pAT  =  0.03 1  s  1  in  the  away  and  toward  architectures, 
respectively.  The  corresponding  amplitudes  Ar  are  0.13°C  and  0.20°C  in  the  away  and  the  toward  architec¬ 
tures,  respectively.  The  ratio  R  =  A(AK)/a(AT)  between  the  jump  in  voltage  (a(AT))  and  temperature 
difference  (A(AT’))  in  the  away  architecture  is  R  =  7.31  mV/  C  . 

Overall,  Figure  3  and  Figure  4  suggest  that  AK(/)  and  AT(/)  evolve  independently,  with  different  rates 
of  increase  such  that  pAV  »  pAT .  The  electric  contribution  affects  only  AK(/).  Thus,  we  conclude  that  the 
electric  and  entropic  contributions  in  the  interaction  of  IR  radiation  with  a  PG  device  are  decoupled.  To  further 
support  this  conclusion,  in  Appendix-1  we  will  compare  the  trends  of  pAV ,  pAT ,  R  =  A(AF)/a(AT)  ,  and 
the  values  of  A(aR)  upon  activating  the  PG  device  with  IR  radiation,  as  done  so  far,  and  conductive  heat 
transfer  from  a  100  Q  resistor  and  a  0.02  A  current.  Only  the  entropic  contribution  is  activated  in  this  case  be¬ 
cause  the  produced  power  of  0.04  W,  corresponding  to  a  temperature  of  «24°C  on  the  PG  device,  is  too  low  to 
produce  significant  blackbody  radiation  to  trigger  the  electric  contribution. 

d)  Behavior  of  AV (/)  and  AT (/)  in  the  30  hours  after  starting  the  illumination 

The  voltage  difference  A V(t)  and  temperature  difference  AT(/)  evolve  in  the  30  hours  after  starting  the 
illumination  of  the  PG  device  as  pictured  in  Figure  5  and  Figure  6.  We  observed  one  possible  behavior  in  Fig¬ 
ure  5,  which  displays  the  data  collected  with  the  toward  architecture.  The  voltage  difference  AV  (/)  in  Figure 
5(a),  exhibits  an  instability  around  the  20th  hour.  No  instability  is  detected  in  A r(/)  in  Figure  5(b),  which  is 
almost  flat.  The  same  instability  affecting  AR(/) ,  however,  is  revealed  in  Tm  (/)  and  Tno_m  (/) ,  displayed  in 
Figure  5(c)  and  Figure  5(d),  respectively.  This  finding  suggests  that  a  strong  correlation  exists  among  AF(/)  , 
Tir  (/)  and  Tno_m(t) ,  which  we  ascribe  to  the  capacitor-type  of  structure  of  the  PG  device,  as  predicted  in  Sec¬ 
tion  3(b). 

We  observed  another  possible  behavior  in  Figure  6,  displaying  data  collected  with  the  away  architecture.  In 
this  case,  AV (/)  ,  shown  in  Figure  6(a),  exhibits  an  instability  around  the  12th  hour  which  is  uncorrelated  with 
the  trends  occurring  in  TJR(t)  and  Tno_m(t).  The  lack  of  correlation  highlights  the  lack  of  coupling  between 
the  electric  and  the  entropic  contributions.  Since  the  trends  of  the  two  temperatures  are  identical,  we  show  only 
Tir  (/)  in  Figure  6(b).  We  omit  to  display  AT(/)  ,  which  is  flat  as  hypothesized  in  Section  3(b).  The  instabili¬ 
ty  of  A f(/)  in  Figure  6(a)  is  non-periodic,  as  predicted  in  Section  3(b).  Thus,  to  find  a  suitable  fitting  func¬ 
tion,  we  recall  the  hypothesis  in  Section  3(b)  where  we  pictured  the  time-dependence  of  AV (/)  in  Figure  6(a) 
as  captured  by  a  hyperbolic  secant  function  with  no  phase  relationship  with  P(i)  and  no  sinusoidal  periodicity. 
Thus,  we  propose: 


AV(t)  =  AV0/f+AV0SCsech 


l~‘c 

H 


(11) 


In  this  expression,  AVoff  and  AVosc  are  the  offset  and  amplitude  of  the  departure  of  A  V(t)  from  the  off¬ 
set,  respectively.  The  positive  or  negative  sign  of  A  Vosc  corresponds  to  a  downward  or  upward  concavity,  re¬ 
spectively,  of  the  instability.  The  critical  time  tc  is  the  instant  in  which  the  maximum  |AKo^|  value  is 
achieved.  Finally,  the  term  H  indicates  the  half  width  at  half  maximum  (HWHM),  or  minimum  (depending  upon 
the  sign  of  AVosc )  of  the  instability.  The  magnitude  of  AVoff  corresponds  to  the  long  term  equilibrium  voltage 
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Fi  gure  5.  Panels  (a),  (b),  (c),  and  (d)  correspond  to  data  eolleeted  from  the  toward  architecture  in  the  whole  time  span  of 
about  30  hours  following  the  start  of  the  illumination  of  the  PG  deviee  with  1R  radiation,  (a)  Voltage  difference  AF(/) .  (b) 

Temperature  difference  AT(/) .  Temperatures  TJR(t)  (c)  and  Tno_m{t)  (d)  of  the  illuminated  and  non-illuminated  faees, 

respectively,  of  the  PG  deviee.  The  zero  of  the  time-seale  eoineides  with  the  start  of  the  illumination  with  the  IR  radiation, 
and  the  unit  for  time  is  the  hour  ([h]). 


A  Vn 


reported  in  Table  3.  The  typical  values  of  A Voff  ,  A Vosc  ,  ic,  and  Fi  are  reported  in  Table  4,  labelled  in 


Figure  6(a),  and  were  obtained  by  placing  the  zero  of  the  time-scale  at  the  start  of  the  illumination.  We  high¬ 


light  that,  since  Tm  (/)  in  Figure  6(b)  does  not  peak  at 


A F(/)  .  The  instability  in  AF(/)  is  evidenced  in  the  AF(/), 


it  does  not  correlate  with  the  behavior  in  time  of 

dAF(/V  f  . 

— graph  in  Figure  6(c),  where  stable 


dt 


and  unstable  fixed  points  [20]  alternate  in  a  complex  fashion  without  periodicity. 

Since  the  time-dependence  is  enclosed  in  a  hyperbolic  secant  function,  we  name  the  instability  in  AF(/)  in 
Equation  (1 1)  and  Figure  6(a)  as  hyperbolic  instability.  We  establish  the  lack  of  correlation  between  A F(/) , 
on  one  hand,  and  Tm  (/)  and  Tno  IR{t)  on  the  other,  as  the  criterion  to  identify  such  instability.  Since  the 
hyperbolic  instability  in  A V[t)  is  absent  in  AT(/) ,  as  expected  from  Section  3(b),  we  relate  the  instability  to 
the  sole  electric  contribution. 

e)  Correlation  between  P{()  and  AF(/)  in  the  30  hours  after  starting  the  illumination 

Here  we  highlight  the  correlations  existing  between  the  power  P(t)  and  the  voltage  difference  AF(/)  to 
further  support  the  choice  of  the  fitting  function  in  Equation  (11)  for  A F(/)  based  on  the  hypothesis  hig¬ 
hlighted  in  Section  3(b). 

1)  The  average  Hp  of  the  power  P(/)  in  Table  2  is  4.0  ±0.5  h,  about  twice  the  value  of  1.82  h  reported 
in  Table  4  for  H  of  AF(/). 

2)  The  instabilities  in  P(t )  and  AF(/)  are  both  small:  the  former  is  Poff  j Posc  « 1 03 ,  while  the  latter  is 


/AKOJC*10/l0 


3)  The  lc  term  in  Table  4  is  found  at  11.65  h,  which  is  reasonably  close  to  the  average  tcP  at  16±2  in 
Table  2. 

These  observations  further  support  our  hypothesis  in  Section  3(b)  relating  the  sinusoidal  instability  of  P{t) 
to  AF(/)  through  on  the  Cu  plate.  In  Appendix-2  we  will  show  that  the  hyperbolic  instability  in 

AF(/),  selected  according  to  the  criterion  established  in  Section  3(d),  occurs  frequently  in  our  observations. 
We  will  provide  further  evidence  of  the  link  between  tc  and  tcP  .  Finally,  we  will  prove  the  lack  of  phase  rela¬ 
tionship  between  F(/)  and  AF(/)  through  lack  of  predictability  of  the  sign  of  the  amplitude  AFovc  of  the 
departure  of  AF(/)  from  equilibrium. 

f)  The  nonlinearity  of  the  hyperbolic  instability  in  AF(/) 

We  found  that  Equation  (11),  used  to  fit  the  hyperbolic  instability  in  AF(/)  in  the  30  hours  after  starting  the 
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AV(t)  [mV] 


(c) 


Figure  6.  Panels  (a),  (b).  and  (c)  correspond  to  data  collected  from  the  away  architecture  in  the  whole  time  span  of  about  30 
hours  following  the  start  of  the  illumination  of  the  PG  device  with  IR  radiation,  (a)  Voltage  difference  AF(/)  exhibiting 

the  hyperbolic  instability  obeying  Equation  (11).  The  parameters  A  Voff ,  A I  ,  tc ,  and  H ,  reported  in  Table  4,  are  labelled, 
(b)  Temperature  Tm  (/)  of  the  face  of  the  PG  device  illuminated  by  the  1R  radiation,  and  exhibiting  no  correlation  with  the 

hyperbolic  instability  in  AT(f)  highlighted  by  tc.  (c)  Graph  of  j  obtained  using  the  parameters  in  Ta¬ 

ble  4.  The  zero  of  the  time-scale  coincides  with  the  start  of  the  illumination  with  the  IR  radiation,  and  the  unit  for  time  is  the 
hour  ([h]). 


Table  4.  Parameters  of  A F(/)  =  AF^  +  AFwcsech|  j  in  Equation  (11)  used  to  fit  the  voltage  difference  AT(/)  in 

Figure  6(a)  collected  from  the  away  architecture  in  the  time  interval  of  about  30  hours  following  the  start  the  illumination  of 
the  PG  device  with  IR  radiation.  The  parameters  AVoff ,  AVmc ,  tc .  and  H  are  labeled  in  Figure  6(a). 

Away  A^[mV]  ATjmV]  t  [h]  H[ h] 

1.037  0.0137  11.65  1.82 


illumination  of  the  PG  device  with  IR  radiation,  can  be  a  solution  of  the  equation: 

,  x  dAV(t)  d3A V(t) 

crAV(t) - 5^  +  c - r^  =  0. 

W  dt  dt 3 


(12) 


Equation  (12)  is  nonlinear  because  of  the  3rd  order  partial  differential  of  A V{t)  with  respect  to  time  t.  The 
equation  resembles  the  Korteweg-de  Vries  (KdV)  equation  [26]  [27]  after  eliminating  the  space-dependence, 
and  flipping  the  space  variable  with  time  /.  In  order  for  AK(/)  in  Equation  (11)  to  be  a  solution  of  Equation 
(12),  we  select  the  o  and  g  coefficients  by  rewriting  Equation  (12)  as: 
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,.dA  V(t)  d}AV  It) 

nAf'O) — — A^=o. 

v  ’  dt  8t 3 


(13) 


By  substituting  Equation  (11)  into  Equation  (13),  we  obtain; 


5sech 


t-t 


v  H 


-tan 


izic 

h 


AKOTCsech 


t-*c 

H 


+  AV 


off 


(14) 


The  time-dependent  coefficient  =  of  Equation  (13)  stands  for  nonlinear  effects  [26],  and 

infers  that  a  complex  dynamics  is  hidden  in  the  hyperbolic  instability  of  A V {t)  [28]  [29],  In  Figure  7  we 
graph  Q(/)  using  the  parameters  reported  in  Table  4  for  the  away  architecture.  Despite  the  complexity  of 
Equation  (14),  Q(/)  appears  surprisingly  simple,  correlates  with  the  time-scale  of  the  hy  perbolic  instability  in 
AV (/)  ,  and  peaks  at  tc.  Furthermore,  we  determined  the  physical  units  of  Q(/)  to  be  l//?2  (I/mV)  (where  h 
indicates  the  hour),  which  suggests  that  Q(/)  could  be  interpreted  as  the  change  in  time  of  the  rate  of 
1/AK(,),  or  simply,  the  “acceleration”  of  the  inverse  voltage.  Such  identification  points  out  the  existence  of 
nonlinear  “forces”  that  continuously  push  and  slow  down  the  changes  in  AF(y)  .  We  plan  on  further  investi¬ 
gating  this  topic. 


4.  Summary  and  Significance 

We  identify  and  decouple  an  electric  and  an  entropic  contribution  to  energy  transfer  from  low  power  infrared 
radiation  to  a  power  generator  device.  The  electric  contribution  is  related  to  the  effects  of  the  electric  E  and 
magnetic  //  fields  in  the  infrared  radiation  and  is  detected  through  the  voltage  produced  by  the  power  genera¬ 
tor  device.  The  entropic  contribution  is  mainly  related  to  the  temperature  difference  between  the  faces  of  the  de¬ 
vice.  Two  observations  enable  us  to  decouple  the  electric  and  entropic  contributions.  First,  the  electric  contribu¬ 
tion  imparts  faster  rates  of  increase  of  the  voltage  difference  AV  (/)  produced  by  the  power  generator  device 
immediately  after  starting  the  illumination  with  infrared  radiation.  Second,  the  electric  contribution  generates  a 
hyperbolic  instability  in  the  30  hours  after  starting  the  illumination.  The  entropic  contribution  changes  slowly 
and  simply  reflects  the  relative  delay  with  which  the  two  faces  of  the  power  generator  device  respond  to  thermal 
perturbations. 

Our  preliminary  studies  suggest  that  it  is  important  to  learn  to  exploit  the  electric  contribution,  because  it  of¬ 
fers  a  variety  of  ways  to  increase  the  voltage  produced  by  the  power  generator  device.  For  example,  black- 
colored  gaffer  tape  on  the  illuminated  face  doubles  the  amount  of  voltage  produced,  as  discussed  in  Appendix- 1, 
while  the  temperature  difference,  thus  the  entropic  contribution,  is  not  sensitive  to  the  presence  of  the  tape.  Our 
findings  are  relevant  for  understanding  the  mechanisms  for  harvesting  1R  radiation,  and  possibly  other  electro¬ 
magnetic  radiations,  through  a  power  generation  device  as  an  alternative  energy  source. 


away  vertical,  0O  =  45° 


Figure  7.  “Acceleration”  of  the  inverse  voltage  Q(/)  for  the  away  architecture  obtained  from  Equation  (14)  using  the  pa¬ 
rameters  in  Table  4.  The  time-dependent  coefficient  Q(/)  peaks  at  tc ,  reported  in  Table  4.  The  unit  for  time  is  the  hour 
([«). 
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The  future  efforts  will  be  devoted  to  overcome  the  limitations  of  the  present  work:  1)  understand  the  relation¬ 
ship  among  infrared  source  power,  surface  charge  density  on  the  power  generator,  and  produced  voltage,  2)  cla¬ 
rify  the  behavior  of  the  power  generator  as  a  capacitor,  and  3)  investigate  the  role  of  infrared  radiation  power  on 
the  results. 
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Appendix-1 

To  highlight  the  effects  of  the  entropic  contribution  and  decouple  it  from  the  electric  contribution,  we  collected 
the  voltage  difference  A V(t)  and  the  temperature  difference  AT(/)  from  the  PG  device  activated  by  con¬ 
ductive  heat  transfer  from  a  1 00  Q  resistor  in  contact  with  the  surface  of  the  device.  A  temperature  of  ~24°C  was 
generated  by  the  PG  device  activated  with  a  B&K  Precision  1665  power  supply  providing  0.02  A  and  3.2  V  to 
the  resistor.  The  achieved  temperature  is  of  the  same  order  of- magnitude  as  that  detected  when  IR  radiation  hits 
the  PG  device,  as  can  be  seen  in  Figure  5(c)  and  Figure  5(d).  We  used  a  PG  device  071 1 1-9L31-04B  by  Cus¬ 
tom  Thermoelectric  Inc.  finished  with  an  AlO  plate  (not  with  a  Cu  plate)  and,  in  selected  cases,  gaffer  tape.  We 
placed  the  sample  holders  on  both  sides  of  the  PG  device  in  the  toward  architecture  to  create  a  capacitor  struc¬ 
ture  with  electrodes  on  both  faces  and  thus  avoid  the  “decay”  in  AK(/)  shown  in  Figure  3(d).  The  experi¬ 
ments  investigating  the  effects  of  conductive  heat  transfer  were  performed  in  an  insulated  sample  compartment 
described  in  Ref.  [30].  The  PG  device  was  horizontally  fixed  on  the  sample  holders  in  all  the  measurements. 
When  the  1R  radiation  was  used  to  compare  the  results  obtained  with  conductive  heat  transfer,  the  PG  device 
was  also  positioned  horizontally  and  at  an  angle  of  incidence  00  =30  with  respect  to  the  IR  radiation.  The  ba¬ 
sic  experimental  parameters  are  summarized  in  Table  1. 

In  the  400  seconds  immediately  following  the  start  the  activation  of  the  PG  device  by  either  conductive  heat 
transfer  or  IR  radiation,  we  measured  the  rate  of  increase  of  A V (t)  {pAV)  and  A T(/)  (pAT)9  and  the  ratio 
R  =  A(AK)/A(Ar)  between  the  jumps  in  voltage  (a(AT))  and  temperature  difference  (a(AT))  in  the 
away  and  toward  architectures.  We  report  the  average  values  of  PAy ,  PAT  >  and  R>  alongside  their  uncertainties, 
in  Figures  Al(a)-(c),  respectively.  Figure  Al(a)  indicates  that  the  average  pAV  values  obtained  by  activating 
the  PG  device  with  IR  radiation  are  always  larger  than  those  obtained  by  activating  the  PG  device  using  conduc¬ 
tive  heat  transfer.  Figure  Al(b)  suggests  that  the  average  pAT  values  have  a  common  value,  independently  of 
the  source  of  activation  of  the  PG  device.  Moreover,  the  average  pAV  and  pAT  values  are  similar  when  the 
PG  device  is  activated  by  conductive  heat  transfer.  In  addition,  the  pAT  values  coincide  with  those  derived  from 
Figure  4(b)  and  Figure  4(d).  This  finding  supports  the  identification  of  the  contribution  to  energy  transfer  with 
slower  rate  of  increase  with  the  entropic  contribution.  Figure  Al(c)  shows  that,  within  the  errors,  the  average 
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Figure  Ah  (a)  Rate  of  increase  of  the  voltage  difference  ( pAV  ).  (b)  Rate  of  increase  of  the  temperature  difference  ( pAT ).  (c) 
Ratio  R  -  A(AV)/a(AT)  between  the  jump  in  voltage  (a(AJ  ))  and  temperature  difference  (A (AT)) .  The  data  were 

obtained  performing  measurements  in  the  away  (//)  and  toward  (T)  architectures  on  a  PG  device  finished  with  an  AlO  plate 
covered  with  gaffer  tape,  described  in  Appendix- 1.  The  results  compare  data  obtained  in  the  400  seconds  immediately  fol¬ 
lowing  the  start  the  activation  of  the  PG  device  with  IR  radiation  (IR)  or  conductive  heat  transfer  (Th)  from  a  100  O  resistor 
The  horizontal  line  at  0.05  s' 1  in  panels  (a)  and  (b)  highlights  the  level  above  which  only  pAV  values  due  to  IR  radiation 
can  be  detected.  The  PG  device  was  horizontally  fixed  on  tire  sample  holders  in  all  the  measurements.  When  the  IR  radiation 
was  used,  the  PG  device  w7as  positioned  at  an  angle  of  incidence  00  =  30°  with  respect  to  the  IR  radiation. 
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value  of  the  ratio  R  is  larger  when  the  PG  device  is  activated  by  1R  radiation  through  the  electric  contribution. 

In  the  400  seconds  immediately  following  the  start  the  activation  of  the  PG  device  by  either  conductive  heat 
transfer  or  IR  radiation,  we  also  measured  the  jump  in  voltage  (a(AF))  obtained  performing  measurements 
on  the  PG  device  described  above  and  finished  with  colored  gaffer  tape.  We  performed  the  measurements  for  the 
away  and  toward  architectures.  Figure  A2(a)  displays  the  magnitude  of  A(AK)  obtained  by  activating  the  PG 
device  with  conductive  heat  transfer  from  a  100  Cl  resistor.  The  observed  trends  are  neither  affected  by  the  color 
nor  by  the  presence  of  the  tape.  On  the  other  hand.  Figure  A2(b)  shows  that  the  magnitude  of  A(AK)  obtained 
by  activating  the  PG  device  with  IR  radiation,  is  slightly  sensitive  to  the  color  of  the  tape,  and  exhibits  a  notice¬ 
able  drop  when  the  tape  is  absent.  Interestingly,  the  black-colored  tape  on  the  illuminated  face  of  the  PG  device 
doubles  the  magnitude  of  A(AK)  compared  to  the  case  without  tape.  Values  of  A(AK)  above  0.5  mV,  repre¬ 
sented  by  the  horizontal  line  in  Figure  A2(b),  can  be  achieved  only  with  tape  present  on  the  illuminated  face  of 
the  PG  device.  Thus,  we  conclude  that  the  tape  profoundly  affects  the  capacitor-type  behavior  of  the  PG  device. 

In  the  time  span  of  about  50  hours  following  the  start  of  the  activation  of  the  PG  device  with  conductive  heat 
transfer  from  a  100  Q  resistor,  the  voltage  difference  AK(/)  ,  temperature  difference  AT (7)  ,  and  temperatures 
Tho({t )  and  Tcold[t)  of  the  activated  and  non-activated  faces,  respectively,  are  shown  in  Figures  A3(a)-(d). 
We  collected  the  particular  set  of  data  chosen  using  the  PG  device  described  above,  finished  with  pink-colored 
gaffer  tape,  and  set-up  in  the  away  architecture.  Nevertheless,  Figure  A3  is  representative  of  all  the  data  col¬ 
lected  in  the  away  and  toward  architectures,  with  tape  of  all  the  available  colors,  and  without  tape.  In  all  these 
cases,  we  observed  A V{t) ,  AT(/),  Thot(t),  and  TcoUJ(t )  to  be  featureless  and  follow  the  same  trends.  The 
trend  of  AF(/)  is  symmetric  to  that  of  AT(/),  7^,  (/),  and  Tcold  (/).  This  behavior  differs  from  that  observed 
for  the  PG  device  activated  by  the  IR  radiation  and  captured  in  Figure  5.  In  this  case,  AV(t)  and  the  tempera¬ 
tures  Tir  (7)  and  Tno_m{t)  of  the  illuminated  and  non-illuminated  faces  are  strongly  correlated,  while 
AT (/)  =  Tir  (/)  -  Tn0_m  (/)  is  flat  and  featureless. 

The  results  in  Figures  A1-A3,  together  with  those  of  Figures  3-5,  further  support  that  the  electric  and  the  en- 
tropic  contribution  to  energy  transfer  from  low  power  IR  radiation  to  the  PG  device  are  decoupled. 

Appendix-2 

Figure  6  shows  the  hy  perbolic  instability  revealed  in  the  voltage  difference  AV  (t}  in  the  time  span  of  about 
30  hours  following  the  start  of  the  illumination  of  the  PG  device  with  IR  radiation  in  the  away  architecture.  We 
identify  such  hy  perbolic  instability  when  no  correlation  is  found  between  AV (/)  and  the  temperatures  T}R  (/) 
and  Tno_m[t )  of  the  illuminated  and  non-illuminated  faces  of  the  PG  device,  respectively,  and  when  AT(V)  is 
flat.  To  demonstrate  that  hyperbolic  instabilities  are  common  in  the  interaction  between  IR  radiation  and  a  PG 
device,  we  collected  a  set  of  measurements  with  the  PG  device  described  in  Appendix- 1.  With  the  AlO  plate  or 


Figure  A2.  Jump  in  voltage  (a(AJ'))  obtained  performing  measurements  in  the  away  (A)  and  toward  (7)  architectures  on 

a  PG  device,  described  in  Appendix- 1,  finished  with  an  AlO  plate  covered  with  colored  gaffer  tape.  The  results  compare  data 
obtained  in  the  400  seconds  immediately  following  the  start  the  activation  of  the  PG  deviee  with  (a)  conductive  heat  transfer 
from  a  100  Q  resistor  or  (b)lR  radiation.  The  PG  deviee  was  horizontally  fixed  on  the  sample  holders  in  all  the  measure¬ 
ments.  When  the  IR  radiation  was  used,  the  PG  deviee  was  positioned  at  an  angle  of  incidence  #0  =  30  with  respect  to  the 
IR  radiation.  The  horizontal  line  at  0.5  mV  represents  the  minimum  value  of  A(AF)  that  can  be  obtained  only  by  plaeing 
tape  on  the  illuminated  face  of  the  PG  deviee. 
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Figure  A3.  Panels  (a),  (b),  (c),  and  (d)  correspond  to  data  collected  from  the  away  architecture  in  the  whole  time  span  of 
about  50  hours  following  the  start  of  the  activation  of  the  PG  device  with  conductive  heat  transfer  from  a  100  Q  resistor,  (a) 
Voltage  difference  A  V  (/) .  (b)  Temperature  difference  A  T  (?) .  Temperatures  7^,(/)  (c)  and  TcoU(t)  (d)  of  the  activated 

and  non-activated  faces,  respectively,  of  the  PG  device.  The  PG  device  was  horizontally  fixed  on  the  sample  holders  in  these 
measurements  and  was  finished  with  a  layer  of  pink-colorcd  gaffer  tape  on  the  AlO  plate.  The  unit  for  time  is  the  hour  (Pi]). 


(b) 


Figure  A4.  (a)  Averages  of  the  critical  time  ( tc )  and  (b)  amplitude  ( A Vwc )  of  the  hyperbolic  instability  detected  in  the  time 

span  of  about  50  hours  following  the  start  of  the  illumination  of  the  PG  device  with  1R  radiation.  The  data  were  collected 
performing  measurements  on  a  PG  device  finished  with  an  AlO  plate,  described  in  Appendix- 1,  and  in  the  away  (A)  and  the 
toward  (7)  architectures.  The  PG  device  was  positioned  horizontally  and  at  an  angle  of  incidence  <90  =305  with  respect  to 
the  IR  radiation,  file  average  of  the  higher  and  lower  values  of  /  found  in  the  away  and  toward  architectures  are  compared 
with  the  average  value  of  tcP  at  1 6  ±  2  h  for  the  power  P[t) .  In  (b),  the  horizontal  line  separates  the  positive  and  nega¬ 
tive  AVose . 


the  layer  of  tape  facing  the  1R  radiation,  the  surface  density  <r(z,t)  refers  to  the  charges  in  the  AlO  plate  or  in 
the  tape.  The  results,  summarized  in  Figure  A4,  report  the  average  values  of  the  critical  time  ( tc )  and  amplitudes 
(AVosc)  of  the  revealed  hyperbolic  instabilities.  We  found  a  correlation  neither  between  the  values  of  tc  in 
Figure  A4(a)  and  the  signs  of  AVosc  in  Figure  A4(b),  nor  between  these  quantities  and  the  away  or  toward 
architectures.  However,  by  observing  Figure  A4(a),  we  found  that  tc  is  typically  located  around  36±lh  , 
which  is  about  twice  the  magnitude  of  tcP  at  16  ±2  h  for  the  sinusoidal  instability  in  P{t)  of  the  IR  source 
reported  in  Table  2.  In  addition,  we  noted  lower  values  for  tc  in  Figure  A4(a)  located  around  20 ±3  h  ,  le.  in 
the  neighborhood  of  tcP .  These  observations  highlight  the  correlations  existing  between  the  instabilities  in 
P(t)  and  A V[t)  as  already  discussed  in  Section  3(e).  Furthermore,  in  Figure  A4(b)  we  noted  randomness  in 
the  sign  of  AVosc ,  which  we  attribute  to  the  loss  of  correlation  between  the  phase  of  P(t)  and  <r(zj ) ,  and 
then  between  a(z,/)  and  AV  (/)  .  We  predicted  the  loss  of  phase  relationship  in  Section  3(b). 

We  summarize  the  findings  in  Figure  A4  as  follows:  the  hyperbolic  instability  in  AV(t)  is  a  common 
phenomenon,  which  is  linked  with  the  sinusoidal  instability  in  the  power  P(t )  of  the  IR  radiation. 


© 


